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As a massive star evolves through multiple stages of nuclear burning on its 
way to becoming a supernova, a complex, differentially rotating structure is set 
up. Angular momentum is transported by a variety of classic instabilities, and 
also by magnetic torques from fields generated by the differential rotation. We 
present the first stellar evolution calculations to follow the evolution of rotating 
massive stars including, at least approximately, all these effects, magnetic and 
non-magnetic, from the zero-age main sequence until the onset of iron-core col- 
lapse. The evolution and action of the magnetic fields is as described by Spruit 
(2002) and a range of uncertain parameters is explored. In general, we find that 
magnetic torques decrease the final rotation rate of the collapsing iron core by 
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about a factor of 30 to 50 when compared with the non-magnetic counterparts. 
Angular momentum in that part of the presupernova star destined to become a 
neutron star is an increasing function of main sequence mass. That is, pulsars 
derived from more massive stars will rotate faster and rotation will play a more 
dominant role in the star's explosion. The final angular momentum of the core 
is determined - to within a factor of two - by the time the star ignites carbon 
burning. For the lighter stars studied, around 15 M©, we predict pulsar periods 
at birth near 15 ms, though a factor of two range is easily tolerated by the un- 
certainties. Several mechanisms for additional braking in a young neutron star, 
especially by fall back, are also explored. 

Subject headings: stars: massive, rotation, magnetic fields - pulsars: rotation 

1. Introduction 

Massive stars are known to be rapid rotators with typical equatorial velocities around 
200kms~^ or more during hydrogen burning (Fukuda 1982). It has long been known that 
this much rotation could have a significant infiuence on the evolution of the star, both on the 
main sequence and during in its later stages of evolution (e.g., Endal & Sophia 1976; Endal & 
Sofia 1978; Heger et al. 2000; Maeder & Mcynct 2001; Hirishi et al. 2004). For example, on the 
main sequence, surface abundance patterns are observed (e.g., Gies & Lambert 1992; Herrero 
1994; Vrancken et al. 2000; Venn 1999) that are most naturally explained by rotationally- 
induced mixing in the stellar interior (e.g., Venn 1999; Heger & Langer 2000; Meynet & 
Maeder 2000). 

Numerical studies of of the later evolutionary stages of rotating massive stars initially 
found that the stellar core would reach critical (Keplerian) rotation by the time carbon 
ignited in the core (e.g., Kippenhahn et al. 1970; Endal & Sophia 1976), especially when 
angular momentum transport is neglected. Thus one might encounter triaxial deformation 
early on, with continued evolution hovering near instability. Modern stellar models that 
include more instabilities capable of transporting angular momentum find that considerable 
angular momentum is lost from the core during hydrogen and helium burning (Endal & Sofia 
1978), especially when the stabilizing effect of composition gradients is reduced (Pinsonneault 
et al. 1989; Heger et al. 2000; Maeder & Mcynct 2001; Hirishi et al. 2004). Such stars now 
evolve to iron core collapse, but still form very rapidly rotating (~ 1 ms) neutron stars. 

However, with the exception of Spruit & Phinncy (1998) and Maeder & Meynet (2004), 
all studies of massive stellar evolution to date have ignored what is probably a major effect. 
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the torques exerted in differentially rotating regions by the magnetic fields that thread them. 
This omission has not been because such torques were thought to be unimportant, but 
because of the complexity and uncertainty of carrying out even one-dimensional calculations 
that included them. This uncertainty, in turn, related to the absence of a credible physical 
theory that would even qualitatively describe the separate evolution of the radial and poloidal 
components of the field, both of which are necessary to calculate the torque. 



Spruit & Phinney (1998) estimated these torques assuming that the poloidal field results 
from differential winding and that the radial field would, due to the action of unspecified 
instabilities, be comparable. Particular attention was paid to the region of large shear that 
separates the core of helium and heavy elements from the very slowly rotating hydrogen 
envelope of a red supergiant star. This assumption, B,. ^ B^. probably overestimated the 
actual Br by orders of magnitude and the rotation rates that Spruit & Phinney (1998) 
estimated for young pulsars were thus too slow (see also Livio & Pringle 1998). 

More recently. Spruit (2002) has provided improved estimates for B^ and B^ based 
upon a dynamo process that takes into account the effect of stable stratifications. Here 
we investigate the effect of this prescription in a stellar evolution code used to study the 
complete evolution of massive stars. Particular attention is given to the rotation rate it 
implies for young neutron stars. In §2, we describe the physical and numerical modeling. In 
§3, the results for 12, 15, 20, 25 and 35 M© stars are presented and, in §4, compared with 
observational data. Typical neutron star rotation rates, at birth, are ~ 15 ms. In §5 we 
discuss other possible ways that such rapidly rotating stars might be slowed during the first 
days of their evolution by magnetic winds and fall back. The surface abundance changes 
due to rotationally-induced mixing are discussed in §6 and in §7 we give our conclusions. 



2. Implementation of Magnetic Braking 

The treatment of magnetic torques given by Spruit (2002) was implemented in a version 
of the imphcit hydrodynamics stellar evolution code (KEPLER, Weaver et al. 1978) that 
already included angular momentum transport and mixing due to a number of non-magnetic 

processes. Sec Hegcr ct al. (2000) for a discussion of these non-magnetic instabilities and 
Spruit (1999, 2002) for a detailed description of the dynamo mechanism. Here we discuss 
only the implementation of and results from that physical model. 

Both angular momentum transport and chemical mixing are calculated by solving the 
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time-dependent diffusion equation (Heger et al. 2000). It is assumed tfiat tfie dynamo adjusts 
and readies equilibrium field strength on a time-scale that is short compared to the evolution 
time-scales of the star, i.e., we neglect the time required to reach the steady state field 
strength and distribution described by Spruit (2002) and only apply their equilibrium values. 
This should be a reasonable approximation in all but the most advanced burning stages that 
occur when the angular momentum distribution has already been frozen in. The validity of 
this assumption is revisited in §7 



2.1. The dynamo process in radiative regions 

Spruit (2002) discusses two limiting cases: stabilization by composition gradients (Spruit's 
Case 0), and stabilization by superadiabatic temperature gradients (Spruit's Case 1). Re- 
alistic stellar models exhibit both kinds of gradients at same time, so the intermediate case 
must also be included (see also Maeder & Mcynet 2004). 

The effective viscosity, Uc, for the radial transport of angular momentum results from 
azimuthal stress due to the field generated by the dynamo, S = BrB^/ifr. Its value is given 
by 

i/e = S/pqQ . (2) 

See also eqs. 34-40 of Spruit (2002). We also implement the "effective diffusivity" , = rje, 
given in his eqs. 41-43. Here fl is the angular velocity, g = dlnfi/dlnr is the shear, r is the 
radius, p is the density, and S is the stress. 

For application outside the regime of "ideal gas with radiation" where the gradient of 
the mean molecular weight determines the stabilization, we replace by a more general 
formulation ("A^omp")) but do not change nomenclature. This way the discussion in Spruit 
(2002) can be followed without further revision. Generally, the Bruntvasala frequency is 

N'^^f-- ^] (3) 
Hp \T, dlnPj ' ^ ' 

where g = —GM[r)/r'^ is the local gravitational acceleration, P is the pressure. Hp = 
— dr/dlnP is the pressure scale height, and the derivatives are total derivatives in the star. 
The compositional contribution to the Bruntvasala frequency is given by 

Hp V dlnP dlnPy' ' ^ ' 

where T is the temperature. The thermodynamic quantities a, 5, and Fi obey their com- 
mon definitions (e.g., Kippenhahn & Weigert 1990). Again, the derivatives are the actual 
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gradients in the star, not thermodynamic derivatives. The thermal contribution to the 
Bruntvasala frequency is then given by 



As stated above, we use the symbol A'^^ to indicate (the square of) the Bruntvasala frequency 
due to the composition gradient for a general equation of state (= A^omp)) ^ot just due to 
changes of mean molecular weight, /i. Appendices A and B give a derivation of the general 
expressions. This formulation is identical with the stability analysis for the different regimes 
(next section) as used in KEPLER, and therefore necessary for consistency in regimes where 
the assumption of an ideal gas is a poor approximation. 



In addition to the radiative regime, I, where both composition gradients and tempera- 
ture gradients are stabilizing, other regimes have to be considered in which "secular" mixing 
processes operate. These are secular as opposed to processes like convection that operate on a 
hydrodynamic time scale. Two of these are semiconvection (II; where a stabilizing composi- 
tion gradient dominates a destabilizing temperature gradient) and thermohalinc convection, 
also called the salt-finger instability (III; where a stabilizing temperature gradient dominates 
a destabilizing composition gradient; e.g., Braun 1997; Heger et al. 2000, their Figure 1 for 
an overview). The final regime of hydrodynamic instability in non-rotating stars is that of 
Ledoux convection (Ledoux 1958; Regime IV; wehere the destabilizing temperature gradi- 
ent is stronger than any stabilization by the composition gradient). In the Ledoux regime, 
effective viscosities and diffusivities are given by mixing length theory (Heger et al. 2000). 
A very approximate treatment of Regimes II and III is given in the next two sections. 



In semiconvective regions (Regime II), the square of the thermal buoyancy frequency 
is negative and thus the limiting Case 1 would be convection. Instead of equation (34) of 
Spruit (2002), we first compute an "dynamo" effective viscosity by 





2.2. The dynamo in semiconvective and thermohaline regions 



2.2.1. The semiconvective case 



(6) 



with i/eo from equation (35) of Spruit (2002) and Qrain = Qo (equations 37-39 of Spruit 2002). 
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Adopting a model for semiconvection in the non-linear regime similar to Spruit (1992), 
we assume the radiative flux is transported by convective motions in the semiconvective 
layers. For the sake of simplicity, we compute an effective turbulent viscosity, z/gc, from 
the velocities according to mixing-length theory (Schwarzschild & Harm 1958) assuming 
Schwarzschild convection (i.e., when the gradient is neglected). 



1/3 

V647rpcprr2 ) ^'^^ 



^ r g<)H,>L Y 



(e.g., as derived from equations 7.6 and 7.7 of Kippenhahn & Weigert 1990) and an assumed 
characteristic length scale of Hp by 

^sc = g-f^P'J^conv, (8) 

(however, see Spruit 1992), and cp is the specific heat capacity at constant pressure. 

Finally, we assume that the actual viscosity is somewhere between that of the case of 
mere semiconvection and that of the dynamo dominated by composition stratification and 
take the geometric mean between the two viscosities, 

i/e = V^rei^sc ■ (9) 

The effective turbulent diffusion coefficient, D^, is computed in an analogous way, Uj-^, i.e., 
Dc = Deofiq), Imin = Qo- Due to the properties of the layers semiconvection model (Spruit 
1992) - negligible diffusivity in typical stellar environments -, we do not compute a geometric 
mean with the semiconvective diffusivity similar to Eq. (9). Instead we just add it to the 
semiconvective diffusion coefficient computed in KEPLER. Typically it is negligible. 



2.2.2. The thermohaline case 

In this regime (III) the compositional buoyancy frequency becomes negative and the 
limiting Case would be Rayleigh- Taylor unstable. We assume the limiting Case 1 of Spruit 
(2002) and compute i^e — J^ei, — Dei, and ^min = analogous to the semiconvective 
case. 

In massive stars, thermohaline convection typically occurs following radiative shell burn- 
ing in the fuel "trace" left behind by an receding central convective burning phase. The 
higher temperature of shell burning can sometimes create heavier ashes than core burning 
and the resulting compositional structure is unstable. The contraction phase after central 
helium burning is the most important domain in massive stars. 
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Generally, the motions associated with the growth of this instabihty ( "salt fingers" ) are 
very slow. Thus the interaction with the dynamo is much more restricted than in the case of 
semiconvection. Therefore we do not employ a interpolation between the dynamo effective 
radial viscosity or diffusivity as is done in §2.2.1. 



2.3. Evaluating the magnetic field strength 

The most interstesting result of the dynamo model is the magnetic field strength. Com- 
bining equations (21), (23), (35) and (36) of Spruit (2002), one has 



i., = r^n^ . (10) 



Using his equation (28), the components of the field are: 

= 167rVVe'g'Qr-2 (11) 
Bi = levrV^i^eg^^^V^ ■ (12) 



3. Results Including Magnetic Torques 

3.1. Implementation 

Based upon the assumptions in §2, the full evolution of stars of 12, 15, 20, 25, and 35 Mq 
of solar metallicity was calculated. For comparison, equivalent models were also calculated 
without magnetic fields for three of these stars. The initial models and input physics were the 
same as in Rauscher et al. (2001), but the total angular momentum of the star was chsen such 
that a typical initial equatorial rotation velocity of ~ 200kms~^ (Fukuda 1982) was reached 
on the zero-age main sequence (ZAMS; see Heger et al. 2000). Time-dependent mixing and 
angular momentum transport were followed as discussed by Heger et al. (2000) with the 
turbulent viscosities and diffusivities from the dynamo model added to those from the model 
for the hydrodynamic instabilities. For the time being, possible interactions between the two 
are neglected (though see Maeder & Meynet 2004). Our standard case (B) implements the 
description as outlined in the previous sections. 



-8- 



3.2. Presupernova Models and Pulseir Rotation Rates 

Tables 1-6 and Fig. 1 summarize our principal results. As also noted by Maeder & 
Meynet (2004), the inclusion of magnetic torques with a magnitude given by Spruit (2002), 
greatly decreases the differential rotation of the star leading to more slowly rotating cores. 

For illustration, we consider in some detail the evolution of the 15 M0 model with 
standard parameter settings. Table 1 shows that, over the course of the entire evolution, the 
total angular momentum enclosed by fiducial spheres of mass coordinate 1.5, 2.5, and 3.5 Mq 
decreases in the magnetic models by a factor of 100 to 200, far more than in the non-magnetic 
comparison model. In both cases most of the angular momentum transport occurs prior to 
carbon ignition (Figs. 2 and 3). The greatest fractional decrease occurs during the transition 
from hydrogen depletion to helium ignition (Table 1, Fig. 2) as the star adjusts to its new 
red giant strTictTirc. During this transition, the central density goes up by a factor of about 
100 (from llgcm^'^ to 1400 gcm^^) with an accompanying increase in differential rotation 
and shear. Though the angular momentum evolves little after carbon burning, modified only 
by local shells of convection, there is a factor of two change after carbon ignition (central 
temperature = 5xlO®K). 

The magnetic torques are sufficiently strong to enforce rigid rotation both on the main 
sequence and in the helium core during helium burning for all masses considered (see also 
Maeder & Meynet 2004). The appreciable braking that occurs during hydrogen and helium 
burning is therefore a consequence of mass loss and evolving stellar structure, especially the 
formation of a slowly rotating hydrogen envelope around a rapidly rotating helium core. 
Typically the more massive stars spend a shorter time during the critical Kelvin-Helmholtz 
contraction between hydrogen depletion and helium ignition and also have a shorter helium 
burning lifetime. Hence, even though they actually have a little less angular momentum at 
the fiducial mass points in Table 1 on the main sequence, the more massive stars have more 
angular momentum in their cores at the end of helium burning. This distinction persists 
and the more massive stars give birth to more rapidly rotating neutron stars. 

Stars that include all the usual non-magnetic mechanisms for rotational mixing and 
angular momentum transport (Heger et al. 2000), but which lack magnetic fields have a 
considerably different evolution and end up with 30 to 50 times more angular momentum in 
that part of their core destined to collapse to a neutron star (Tables 1 and 3, Figs. 2 and 
3). Most notable is the lack of appreciable braking of the helium core by the outer layers in 
the period between hydrogen core depletion and helium ignition. The post-carbon-ignition 
braking that accounted for an additional factor of two in the magnetic case is also much 
weaker. In fact, the inner 1.5 Mq of the 15 Mq model has only about five times less angular 
momentum at core collapse than implied by strict conservation from the (rigidly rotating) 
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main sequence onwards. As noted previously many times, such a large amount of angular 
momentum would have important consequences for the supernova explosion mechanism. 
Rigidly rotating neutron stars cannot have a period shorter than about 1 ms. 

Table 2 gives the magnetic field components that exist in different parts of the star 
during different evolutionary phases. The magnetic field is highly variable from location 
to location in the star and the values given are representative, but not accurate to better 
than a factor of a few (Fig. 1), especially during the late stages. As expected, the toroidal 
field, which comes from differential winding is orders of magnitude larger than the radial 
field generated by the Taylcr instability. The product B^Br scales very approximately as 
the central density with the ratio Br/B^ ~ 10~^ — 10~^. Though the formulae used here are 
expected to break down during core collapse, a simple extrapolation to 10^^ g cm~^ for the 
central density of a neutron star suggests ~ 10^^ and Br ~ 10^°. 

Table 3 gives the expected pulsar rotation rates based upon a set of standard assump- 
tions. Following Lattimer & Prakash (2001), it is assumed that all neutron stars (with 
reahstic masses) have a radius of 12 km and a moment of inertia 0.35 Mgrav-R^- Adding the 
binding energy to the gravitational mass, assumed in Table 3 to be 1.4 Mq gives the baryonic 
mass of the core that collapsed to the neutron star, 1.7 M©. Taking the angular momen- 
tum inside 1.7 Mq from the presupernova model, assuming conservation during collapse, and 
solving for the period then gives the values in the table. 

3.3. Sensitivity To Uncertain Parameters 

It is expected that our results will be sensitive to several uncertain factors including the 
effect of composition gradients, the efficiency of the dynamo in generating magnetic field, and 
the initial angular momentum of the star. To some extent, these uncertainties arc ameliorated 
by the strong sensitivity of the expected magnetic torques to rotation speed and differential 
shear (eq. 34-36 of Spruit 2002). Composition gradients have very significant infiuence on the 
resulting angular momentum transport by inhibiting the action of the dynamo. In particular, 
composition interfaces, as resulting from the nuclear burning processes in combination with 
convection, usually also show significant shear, but if the dynamo action is suppressed by 
a steep change in mean molecular weight, angular momentum can become "trapped" in a 
fashion similar to that observed by Heger et al. (2000). To illustrate the uncertainties of 
the dynamo model, we studied several test cases where the stabilizing effect of composition 
gradients is varied by multiplying iV^ and by 1/10. We also explored multiplying the 
overall coefficient of the torque in eq. (1) by varying B^j^Br by a factor of ten and multiplied 
the initial rotational rate by 0.5 and 1.5 (Table 3). 
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in three cases, the 20 and 25 M© with high rotation rate, and the "normal" 35 Mq model, 
the star lost its hydrogen-rich envelope due to mass loss either during (25 and 35 Mq) or just 
at the end (20 Mq) of it evolution. The KEPLER code was not able to smoothly run these 
models through the period where they lost their last solar mass of envelope because violent 
pulsations were encountered. The last 0.5 to 1.0 Mq was thus removed abruptly and the 
calculation continued for the resulting Wolf-Rayet star. However, it is clear that removing 
the envelope removes an appreciable braking torque on the helium core so it is reasonable 
that the remnants of such stars rotate more rapidly. 



3.4. Variable Neutron Steir Masses 

Since the iron core masses, silicon shell masses, and density profiles all differ appreciably 
for presupernova stars in the 12 to 35 Mq range, it is not realistic to assume that they all make 
neutron stars of the same mass. Ideally, one would extract realistic masses from a reliable, 
credible model for the explosion, but, unfortunately, such models have not reached the point 
of making accurate predictions. Still there are some general features of existing models 
that provide some guidance. As noted by Weaver et al. (1978), the remnant mass cannot 
typically be less than the mass of the neutronized iron core. Otherwise one overproduces rare 
neutron-rich isotopes of the iron group. More important to the present discussion, successful 
explosions, when they occur in modern codes, typically give "mass cuts" in the vicinity 
of large entropy jumps. The reason is that a jump in the entropy corresponds to sudden 
decrease in the density. This in turn implies a sudden fall off in accretion as the explosion 
develops. The lower density material is also more loosely bound by gravity. The largest 
jump in the vicinity of the iron core is typically at the base of the oxygen burning shell and 
that has sometimes been used to estimate neutron star masses (Timmes et al. 1996). Recent 
studies by Thomas Janka (private communication) suggest a numerical criterion, that the 
mass cut occurs where the entropy per baryon equals AkB- In Table 4, instead of assuming 
a constant baryonic mass of 1.7 Mq for the neutron star progenitor (as in Table 3), we take 
the value where the specific entropy is 4A;B/baryon. 

As expected this prescription gives larger baryonic masses for the remnant ( "Baryon" in 
the table). The fraction of this mass that is carried away by neutrinos is given by Lattimer 
& Prakash (2001), 

f - (13) 

where P — GMgra.v/Rc^- When this is subtracted, one obtains the gravitational masses, 
Mgrav in Table 4. Using these, assuming once more a moment of inertia / = 0.35Mgrav-R^ 
and a radius of 12 km, and conserving angular momentum in the collapse one obtains the 
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period. Because the angular momentum per unit mass increases as one goes out in the star, 
using a larger mass for the remnant increases its rotation rate. 

The resulting values are similar to those in Table 3 for the same "standard" parameter 
settings, but show even more clearly the tendency of larger stars to make more rapidly 
rotating pulsars. They also show an additional prediction of the model - that more massive 
pulsars will rotate more rapidly. The numbers in Table 4 have not been corrected for the 
fact that the typical neutrino, since it last interacts at the edge of the neutron star, carries 
away more angular momentum than the average for its equivalent mass. Thus the periods 
in Table 4 can probably be multiplied by an additional factor of ~ 1.2 (Janka 2004). 

4. Comparison With Observed Pulsar Rotation Rates 

Table 5 gives the measured and estimated rotation rates and angular momenta of several 
young pulsars at birth (Muslimov & Page 1996; Marshall el al. 1998; Kaspi et al. 1994). 
To estimate the angular momenta in the table, we have assumed, following Lattimer & 
Prakash (2001), that the moment of inertia for ordinary neutron stars (not quark stars) is 
/ ~ 0.35Mi?^ pa 1.4x 10"^^ g cm^. Here a constant fiducial gravitational mass of 1.4 is 
assumed. 

Given the uncertainties in both our model and the extrapolation of observed pulsar 
properties, the agreement, at least for the 12 - ISM© models, is quite encouraging. This 
true all the more so when one reahzes that the periods in Tables 3 and 4 should be multiplied 
by a factor of approximately 1.2 — 1.3 to account for the angular momentum carried away by 
neutrinos when the neutron star forms (Janka 2004). Then our standard, numerically typical 
supernova, 15 M©, would produce a neutron star with angular momentum 5.8x10^^ ergs 
("PreSN" entry in Table 1 divided by 1.3). It should also be noted that the most common 
core collapse events - by number - will be in the 12 to 20 Mq range because of the declining 
initial mass function. Moreover, stars above around 25 - 30 Mq may end up making black 
holes (Fryer 1999). 

What then is the observational situation? Marshall el al. (1998) have discovered a pulsar 
in the Crab-like supernova remnant, N157B, in the Large Magellanic Cloud, with a rotation 
rate of 16 ms. This is probably an upper bound to its rotation rate at birth. Glendenning 
(1996) gives an estimated initial rotation rate for the Crab pulsar itself of 19 ms. Estimating 
the initial rotation rate of the much larger sample of pulsars, known to be rotating now 
much more slowly, is fraught with uncertainty. Nevertheless, PSR 0531+21 and PSR 1509- 
58 are also estimated to have been born with ~ 20 ms periods Mushmov & Page (1996) 
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(though 0540-69 is estimated to have Pq 39 ms). At the other extreme, Pavlov et al. 
(2002) find for PKS1209-51/52 a rotation period of 424 ms and argue that the initial rate 
was not much greater. Similarly, Kramer et al. (2003) argue that the initial rotation rate of 
PSR J0538+2817 was 139 ms. 



5. Further Braking Immediately After the Explosion 

5.1. r-Mode Instability 

For some time it was thought that gravitational radiation induced by the "r-mode 
instability" would rapidly brake young neutron stars. However, Arras et al. (2003) found 
that the r-mode waves saturate at amplitudes lower than obtained in previous numerical 
calculations (Lindblom et al. 2001) that assumed an unrealistically large driving force. Much 
less gravitational radiation occurs and the braking time for a rapidly rotating neutron star 
becomes millennia rather than hours. Arras et al. calculate that 

Tspin down = 2000 yr i-^r' {——y\ (14) 

with CKe, the maximum amplitude of the instability (for which Arras et al. estimate 0.1 as 
an upper bound). A 3ms pulsar will thus take about 2000 years to substantially brake and 
the more slowly rotating neutron stars in Tables 3 and 4 will take even longer. 

5.2. Neutrino-Powered Magnetic Stellar Winds 

Woosley & Hcger (2004) and Thompson et al. (2004) have discussed the possibility 
that a young, rapidly rotating neutron star will be braked by a magnetic stellar wind. This 
wind, powered by the neutrino emission of the supernova, may also be an important site for 
r-process nucleosynthesis (Thompson 2003). 

During the first 10 s of its fife a neutron star emits about 20 % of its mass as neutrinos. 
This powerful fiux of neutrinos passing through the proto-neutron star atmosphere drives 
mass loss and the neutron star loses ~ 0.01 M0 (Duncan et al. 1986; Qian & Woosley 1996; 
Thompson et al. 2001). Should magnetic fields enforce corotation on this wind out to a radius 
of 10 stellar radii, appreciable angular momentum could be lost from a 1.4 M© neutron star 
(since j ~ r'^cu). The issue is thus one of magnetic field strength and its radial variation. 

For a neutron star of mass M and radius, Rq, in units of 10 km, the mass loss rate is 
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approximately (more accurate expressions are given by Qian & Woosley) 

M ^ 10-=^ Mo s-^ ( -^^^X" Re'/' ( V ■ (15) 
° VlO ergs/ V ^ / 

The field will cause corotation of this wind out to a radius (Mestcl & Spruit 1987) where 

P (^wind + ^'r') ^ ^. (16) 

Calculations that include the effect of rotation on the neutrino-powered wind have yet to be 
done, but one can estimate the radial density variation from the one-dimensional models of 
Qian & Woosley (1996). For a mass loss rate of 10^^ Mq s'^ they find a density and radial 
speed at 100 km of lO^gcm"^ and 1000 km s~^ respectively. At this time, the proto- neutron 
star radius is 30 km and these conditions persist for ~ 1 s. Later, they find, for a mass 
loss rate of IO^^MqS"^ and a radius of 10 km, a density at 100 km of ~ lO^gcm"^ and 
velocity 2000 km s^^. This lasts ~ 10 s. Unless oo is quite low, Vwind is not critical. For 
u ~ 1000 rads^^ the field required to hold lO^gcm"^ in corotation at 100 km is ~ 3xlO^^G; 
for lO^gcm"^ it is 10^^ G. 

To appreciably brake such a rapidly rotating neutron star, assuming B ~ r^^, thus 
requires ordered surface fields of ~ 10^^ - 10^^ G. Similar conclusions have been reached 
independently by Thompson (2003). Such fields are characteristic of magnetars, but probably 
not of ordinary neutron stars. On the other hand if the rotation rate were already slow, 
u; ~ lOOrads"^ (P ~ 60 ms), even a moderate field of 10^^ G could have an appreciable 
effect. It should also be kept in mind that the field strength of a neutron star when it is 1 - 
10 s old could be very different than thousands of years later when most measurements have 
been made. 



5.3. Fallback and the Propeller Mechanism 

After the first 1000 s, the rate of accretion from fall back is given (MacFadyen et al. 
2001) by 

10-^ ^5 M© s-\ (17) 
with the time in units of 10^ s. For the mass loss rate in units of lO^^gs"^, we obtain 

M26 ~ 2t5-'/=^gs-^ (18) 

For a dipole field with magnetic moment fi^Q = BuRq"^, with B12 the surface field in units 
of 10^^ G, the infalling matter will be halted by the field at the Alfven radius (Alpar 2001), 

TA = 6.8 1^%' M^^" km. (19) 
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At that radius matter can be rotationally ejected, provided the angular velocity there corre- 
sponding to co-rotation exceeds the Keplerian orbital speed. The ejected matter carries away 
angular momentum and brakes the neutron star. This is the propeller mechanism (Illarionov 
& Sunyaev 1975; Chevalier 1989; Lin et al. 1991; Alpar 2001, but see the critical discussion 
in Rapport et al. 2004). 

Obviously ?x must exceed the neutron star radius (10 km here) if the field is to have 
any effect. The above equations thus require a strong field, B > 10^^ G, and accretion rates 
characteristic of an age of at least one day. Additionally there is a critical accretion rate, 
for a given field strength and rotation rate, above which the co-rotation speed at the Alfven 
radius will be slower than the Keplerian orbit speed. In this case the matter will accrete 
rather than be ejected. Magnetic braking will thus be inefficient until o;^ > CM, or 

M26 < 5.7x10"^ c^g/^ /x^o, (20) 

where cu^ is the angular velocity in thousands of radians per second (a;3 = 1 implies a period 
of 27r ms). This turns out to be a very restrictive condition. For a given /i^o and cus, Eqs. (18) 
and (20) give a time, t5,min, when braking can begin. The torque on the neutron star from 
that point on will be 

loj = i_iyr^ = 10^°^, (21) 

with /, the moment of inertia of the neutron star, approximately 10"^^ (Lattimer & Prakash 
2001). The integrated deceleration will be 

Au = 250i,%h-'^,. (22) 

Putting it all together, a 1.4 M© neutron star can be braked to a much slower speed 
{Auj ~ uj) by fall back if fi^o > 78, 43, or 25 for initial periods of 6, 21, or 60 ms, respectively. 
If the surface field strength - for a dipole configuration is less than 2xlO^^G, braking by the 
propeller mechanism will be negligible in most interesting situations. 

However, we have so far ignored all non-magnetic forces save gravity and centrifugal 
force. A neutron star of age less than one day is in a very special situation. The accretion 
rate given by Eq. (17) is vastly super-Eddington. This means that matter at the Alfven radius 
will be braked by radiation as well as centrifugal force and the likelihood of its ejection is 
greater. Fryer et al. (1996) have considered neutron star accretion at the rates relevant here 
and find that neutrinos released very near the neutron star actually drive an explosion of the 
accreting matter. Just how this would all play out in a multi-dimensional calculation that 
includes magnetic braking, rotation, and a declining accretion rate has yet to be determined, 
but is worth some thought. 
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If the accreting matter is actually expelled under 'propeller' conditions, one might expect 
that the mechanism also inhibits the accretion so the process might be self-limiting. This 
traditional view of the propeller mechanism may be misleading however. There is a range 
in conditions where accretion accompanied by spindown is expected to occur (Sunyaev & 
Shakura 1977; Spruit & Taam 1993; for a recent discussion see Rapport et al. 2004), as 
observed in the X-ray pulsars. If fall back is the way most neutron stars are slowed, one 
might expect a correlation of pulsar period with the amount of mass that falls back and that 
might increase with progenitor mass. 

More massive stars may experience more fall back (Wooslcy & Weaver 1995) and make 
slowly rotating neutron stars. But, on the other hand, as Tables 3 and 4 show, neutron stars 
derived from more massive stars are born rotating more rapidly. There may be a mass around 
15 Mq or so where the two effects combine to give the slowest rotation rate. Interestingly, 
the Crab pulsar was probably once a star of ~ 10 M© (Nomoto et al. 1982) and may thus 
have experienced very little fall back. The matter that is ejected by the propeller mechanism 
in more massive stars could contribute appreciably to the explosion of the supernova and 
especially its mixing. 

5.4. Effect of a Primordial Field 

As noted above, the results agree with the initial rotation periods found in Crab-like 
pulsars, but not with the long periods inferred for pulsars like PSR J0538-I-2817. Within the 
framework of our present theory, slowly rotating pulsars must have formed from stars with 
different or additional angular momentum transport mechanisms. Another possibility is that 
their main sequence progenitors started with a qualitatively different magnetic field config- 
uration. The dynamo process envisaged here assumes that the initial field of the star was 
sufficiently weak, such that winding-up by differential rotation produces the predominantly 
toroidal field in which the dynamo process operates. Any weak large scale field initially 
present is then expelled by turbulent diffusion associated with the dynamo. 

Stars with strong initial fields, such as seen in the magnetic A-stars, cannot have followed 
this path. In these stars, the magnetic field is likely to have eliminated any initial differential 
rotation of the star on a short time scale (Spruit 1999), after which it relaxed to the observed 
stable configurations (Braithwaite & Spruit 2004). In these configurations the radial and 
azimuthal field components are of similar magnitude. To the extent that such magnetic 
fields also exist in massive MS stars (where they are much harder to detect, see however 
Donati et al. 2002), they are likely to have led to a stronger magnetic coupling between core 
and envelope, and hence to more slowly rotating pre-SN cores. 
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6. Surface Abundances 

The abundances of certain key isotopes and elements on the surfaces of main sequence 
stars and evolved supergiants are known to be diagnostics of rotationally-induced mixing 
(Heger & Langer 2000; Maeder & Meynet 2000). In particular, on the main sequence, one 
expects enhancements of ^He, ^^C, ^"'N, ^^O, and ^^Na and extra depletion of ^^C, ^^N, ^^O, 
18 19 g ^gg^ these predictions and their sensitivity to assumed magnetic 
torques in the present models, we calculated three additional versions of our 15 M© model in 
which a large nuclear reaction network was carried in each zone. The adopted nuclear physics 
and network were essentially the same as in Rauscher et al. (2001). One model assumed no 
rotation; a second assumed rotation without magnetic torques; and a third, which included 
both rotation and magnetic torques, was our standard model in Tables 1-4. The initial 
rotation rate assumed in the magnetic and non-magnetic models was the same. 

The results (Table 6) show the same sort of rotationally-induced enhancements and 
depletions on the main sequence (defined by half-hydrogen depletion at the center) as in 
Heger & Langer (2000). Helium is up by about half a percent and nitrogen is increased by 
two. Carbon is depleted by about 25 % and ^^N is depleted by a factor of two. The isotopes 
^^C and ^^O are enhanced by factors of three and 1.6 respectively. Be and B are depleted. 

Most important to the current discussion, we see no dramatic differences between the 
surface abundances calculated with and without magnetic torques. Maeder & Meynet (2004) 
suggested that such differences might be a discriminant, perhaps ruling out torques of the 
magnitude suggested by Spruit (2002). 

7. Conclusions and Discussion 

Given the uncertainty that must accompany any first-principles estimate of magnetic 
field generation in the interior of stars, the rotation rates, we derive here compare quite 
favorably with those inferred for a number of the more rapidly rotating pulsars. Without 
invoking any additional braking during or after the supernova explosion, the most common 
pulsars, which come from stars of 10 - 15 M©, will have rotation rates around 10 - 15 ms 
(Table 3). Reasonable variation in uncertain parameters could easily increase this to 20 ms. 
Given that the torques in the formalism of Spruit (2002) depend upon the fourth power 
of the shear between layers, this answer is robust to small changes in the overall coupling 
efficiency and inhibition factors. If we accept the premise that only the fastest solitary 
pulsars represent their true birth properties and the others have been slowed by fall back 
or other processes occurring during the first few years of evolution, our predictions are in 
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agreement with observations. 

If correct, this conclusion will have important implications, not only for the evolution 
of pulsars, but for the supernova explosion mechanism, gravitational wave generation, and 
for gamma-ray bursts (GRBs). For rotation rates more rapid than 5 ms, centrifugal forces 
are an important ingredient in modeling the collapse and the kinetic energy available from 
rotation is more than the 10^^ erg customarily attributed to supernovae. At 10 ms, the energy 
is nearly negligible. 

Given recent observational indications that supernovae accompany some if not all GRBs 
of the "long-soft" variety (e.g., Woosley & Bloom 2005), rotation is apparently important 
in at least some explosions. Otherwise no relativistic jet would form. Indeed, the only 
possibilities are a pulsar with very strong magnetic field and rotation rate ~ 1ms (e.g.. 
Wheeler, Yi, Hoflich, & Wang 2000) and a rapidly rotating black hole with an accretion disk 
(Woosley 1993; MacFadyen & Woosley 1999). Both require angular momenta considerably 
in excess of any model is Table 3 except those that did not include magnetic torques. 

We shall address this subject in a subsequent paper, but for now note a possible im- 
portant symmetry braking condition - the presence of a red supergiant envelope. The most 
common supernovae - Type Hp - will result from stars that spent an extended evolutionary 
period with a rapidly rotating helium core inside a nearly stationary hydrogen envelope. 
GRBs on the other hand come from Type Ic supernovae, the explosions of stars that lost 
their envelopes either to winds or binary companions. Whether a comprehensive model can 
be developed within the current framework that accommodates both a slow rotation rate for 
pulsars at birth and a rapid one for GRBs remains to be seen but we are hopeful. 

As noted in §2, the magnetic torques have been computed assuming that the dynamo 
process is always close to steady state. The time to reach this steady state scales as some 
small multiple of the time for an Alfven wave to travel around the star along the toroidal 
field. This time is very short throughout most of the evolution, but in the latest stages it 
eventually becomes longer than the time scale on which the moment of inertia of the core 
changes. After this, the dynamo effectively stops tracking the changing conditions and the 
field is frozen in (its components varying as 1/r^, for homologous contraction). With the 
data in Tables 1 and 2, we can estimate the effect this would have on the angular momentum 
loss from the core. The Alfven crossing time first exceeds the evolution time around the end 
of carbon burning. Evolving the field under frozen-in conditions from this time on then 
turns out to produce field strengths that do not differ greatly (less than a factor of 4) from 
those shown in Table 2. Since the angular momentum of the core does not change by more 
than 5 % after carbon burning, the effect of using frozen field conditions instead of a steady 
dynamo is therefore small. 
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A. Thermodynamic derivatives 



Instead of the formulations 



1 



1-1 

r2 



Vad = 



7^ 



for an ideal gas with radiation {TZ is the gas constant; Kippenhahn & Weigert 1990), the 
general expressions should be used for more general equations of state, as to, e.g., include 
the effect of degeneracy which is important in the late stages of stellar evolution of massive 
stars and in low mass stars. The thermodynamic derivatives at constant entropy 

'dlnP^ 



Tl = Tad 



1 - - = Vad 
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ad 
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- 1 r2 

(Kippenhahn & Weigert 1990) can be derived from ds = dq/T, dq 
du — P/ dp, and the total derivative on the specific internal energy 



du 

assuming adiabatic changes, ds = 0: 

dq 1 



= ds 



T T 



fdu\ _P_ 



dp + 
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(Kippenhahn & Weigert 1990). This can be transformed into 
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and we readily obtain 
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(A9) 



as a function of the partial derivatives of u with respect to T and p. Using the equation of 
state for bubble without mixing or composition exchange with its surrounding. 
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where 



a 



dlnp 
dlnP 



dlnp 
dlnT 



T fdP 



p_ (dP\ ■ 

p \dp)^_ 



'P r /p, 

and using the thermodynamic relation (e.g., Kippenhahn & Weigert 1990) 

'dP\ P5 



dT 



Ta 



to transform the last term, we can write Eq. (A7) as 



0^ds=i- 



dp A 



P 
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and obtain 

dlnP^ 
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as a function of the derivative of u and P with respect to T and p only. Using Eq. 
can solve for r2 employing the relation for V^: 
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B. Compositional Brunt vasala frequency 

We derive the compositional fraction of the Bruntvasala frequency from the following 
consideration: we compare a change of density if we move a fluid element to a new location 
and allow adjusting it pressure and temperature to its surroundings. When moving to the 
new location, P will be different by dP and T be different by dT. The density change we 
then obtain from the equation of state, Eq. (AlO), for a displacement with such a change in 
pressure thus is 

^^"^^ ^a-6(^\ (Bl) 



dlnPAlement VdlnP/, 

The buoyancy is obtained by comparing with a density change in the star over the same 
change in pressure, the remainder than is due to compositional changes. 

dln(Ap) \ ^ _ dlnT _ d hi p 
dlnP " dlnP dlnP ^ ' 

/ comp 



- 21 - 



For a compositionally stably stratified medium this quantity is negative since pressure de- 
creases outward. 

As an example let us consider an ideal gas with radiation. For the surrounding, the 
equation of state is, as a function of T, P, and mean molecular weight, /i, given by 

dp dP j-dT dp 

— = a— — d—— + (f — (B3 

with 

Substituting this in for the last term of Eq. (B2), we obtain 

dln(Ap)\ dlnu ^ . ^. 

/ comp 

and the usual expression for buoyancy in an ideal gas with radiation is recovered. 

For a displaced blob of material the square of the oscillation "angular" frequency 
{u — 1/period), the Bruntvasala frequency, can be obtained by multiplication with g/Hp, 
analogous to the harmonic oscillator, here for the case of isothermal displacements that only 
consider stabilization by composition gradients: 

?\r2 9 f rdlnT dlnp\ 
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Table 1: Evolution of Angular Momentum at Fiducial Mass Coordinates for ISM© star 



.magnetic star non-magnetic star. 





J(1.5) 


J(2.5) 


J(3.5) 


J(1.5) 


J(2.5) 


J(3.5) 


ZAMS 


1.75x10™ 


4.20x10™ 


7.62x10™ 


2.30x10™ 


5.53x10™ 


1.00x10^1 


H-burn" 


1.31x10™ 


3.19x10™ 


5.83x10™ 


1.51x10™ 


3.68x10™ 


6.72x10™ 


H-dep^ 


5.02x10^9 


1.26x10™ 


2.37x10™ 


1.36x10™ 


3.41x10™ 


6.37x10™ 


He-ign'^ 


4.25x10^8 


1.21x10^9 


2.57x10^9 


1.16x10™ 


2.98x10™ 


4.87x10™ 


He-burn"' 


2.85x10^^ 


7.84x10^^ 


1.83x10^9 


7.06x10^^ 


1.85x10™ 


3.86x10™ 


He-dep"^ 


2.23x10^^ 


5.95x10^^ 


1.21x10^9 


4.72x10^9 


1.26x10™ 


2.52x10™ 


C-ign-'^ 


1.88x10^8 


5.52x10^^ 


1.12x10^9 


4.69x10^'' 


1.26x10™ 


2.46x10™ 


C-dep5 


8.00x10^^ 


3.26x10^^ 


9.08x10^^ 


4.06x10^'' 


1.25x10™ 


2.24x10™ 


0-dep'* 


7.85x10^^ 


3.19x10^^ 


8.43x10^^ 


3.94 xlO^'' 


1.20x10™ 


1.99x10™ 


Si-dep* 


7.76x10^^ 


3.05x10^^ 


7.23x10^^ 


3.75x10^9 


1.16x10™ 


1.95x10™ 


PreSN^' 


7.55x10^^ 


2.59x10^^ 


7.31x10^8 


3.59x10^9 


1.09x10™ 


1.94x10™ 



NOTE: "40 % central hydrogen mass fraction; ''1 % hydrogen left in the core; '^1 % helium 
burnt; "^50% central helium mass fraction; ^1 % helium left in the core; -^central temperature 
of 5xlO^K; ^central temperature of 1.2xlO^K; ''central oxygen mass fraction drops below 
5%; 'central Si mass fraction drops below 10~^; •'infall velocity reaches 1000 km s~^. 
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Table 2: Approximate Magnetic Field and Angular Velocity Evolution in a a 15 M0 Star 



evolution stage 


(Mq) 




(10^ cm) 


(G) 


(G) 


n 

(rads-1) 


MS" 


5.0 


1.9 


90 


2x10^ 


0.5 


4xl0~^ 


TAMS*" 


3.5 


2.7 


67 


3x10^ 


1 


2x10-^ 


He ignition*^ . . 


3.5 


0.90 


52 


5x10^ 


5 


4x10-^ 


He ignition'^ . . 


1.5 


470 


9.8 


3x10^ 


20 


4x10-5 


He depletion"^ 


3.5 


150 


14 


2x10^ 


3.5 


3x10-5 


C ignition^ . . . 


1.5 


1x10^ 


3.5 


6x10^ 


250 


2x10-^ 


C depletion-^ . 


1.2 


3x10^ 


0.80 


3x10^ 


5x10-'^ 


1x10-3 


depletion^ . 


1.5 


4x10^ 


0.73 


2x10^ 


2x10^ 


1x10-3 


Si depletion'* . 


1.5 


2x10^ 


0.44 


5x10^ 


5x10=^ 


3x10-3 


pre-SN^ 


1.3 


5x10^ 


0.12 


5x10*^ 


1x10^ 


5x10-2 


ovohition stage 


(M.) 


Pc 

Vo / 


Tc 

(10" K) 








MS 


5 n 


5 fi 


035 


2.0x10^^ 






TAMS 


3.5 


11 


0.045 


6.4x10^3 






He ignition . . . 


3.5 


1400 


0.159 


6.0x10^3 






He ignition . . . 


1.5 


1400 


0.159 


6.0x10^^ 






He depletion . 


3.5 


2700 


0.255 


1.4x10^2 






C ignition .... 


1.5 


3.8x10^ 


0.50 


2.4x10" 






C depletion . . 


1.2 


7.0x10^ 


1.20 


3.4x10^ 






depletion . . 


1.5 


1.0x10^ 


2.20 


1.1x10^ 






Si depletion . . 


1.5 


4.8x10^ 


3.76 


8.3x10^ 






pre-SN 


1.3 


8.7x10^^ 


6.84 


0.5 







NOTE: 



"40 % central hydrogen mass fraction; 

^1 % hydrogen left in the core; 

^1 % heUum burnt; 

'^1 % hehum left in the core; 

^central temperature of 5xlO®K; 

■^central temperature of 1.2xlO^K; 

^central oxygen mass fraction drops below 5 %; 

'^central Si mass fraction drops below 10"^; 

*infall velocity reaches 1000 km s"^. 
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Table 3: Pulsar Rotation Rate Dependence on Dynamo Model Parameters" 

— Nf^ — — — .B^Br - -VtzAMS- 

0.1 10 0.1 10 0.1 10 0.5 1.5 B=0 
, period (ms) 



initial 


"std." 


mass 




12 M© 


9.9 


15 Mq 


11 


20 M0 


6.9 


25 Mq 


6.8 


35 M . 


1.1'' 



24 4.4 12 10 5.7 21 9.8 10 0.20 
14 3.2 8.4 6.4 3.3 11 7.2 6.5*^ 0.21 
13 3.1 7.3 4.9 2.6 13 7.1 4.3'' 0.22 



"All numbers here can be multiplied by 1.2 to 1.3 to account for the angular momentum 
carried away by neutrinos 

^Became a Wolf-Rayet star during hehum burning 



Table 4: Pulsar Rotation Rate With Variable Remnant Mass" 



Mass 


Baryon^ 


GravitationaP 


J(Mbary) 


BE 


Period*^ 




(Mo) 


(Mo) 


(10^^ ergs) 


(1053 gj,g^ 


(ms) 


12 M© 


1.38 


1.26 


5.2 


2.3 


15 


15 M© 


1.47 


1.33 


7.5 


2.5 


11 


20 M© 


1.71 


1.52 


14 


3.4 


7.0 


25 M© 


1.88 


1.66 


17 


4.1 


6.3 


35 M© ^ 


2.30 


1.97 


41 


6.0 


3.0 



"Assuming a constant radius of 12 km and a moment of inertia 0.35Mi?^ (Lattimer & Prakash 
2001) 

''Mass before collapse where specific entropy is 4/cB/baryon 
■^Mass corrected for neutrino losses 

"^Not corrected for angular momentum carried away by neutrinos 
^ Becaame a Wolf-Rayet star during helium burning 



Table 5: Periods and Angular Momentum Estimates for Observed Young Pulsars 

current initial Jo 



PSR J0537-6910 (N157B, LMC) 

PSR B0531+21 (crab) 

PSR B0540-69 (LMC) 

PSRB1509-58 



16 ~10 8.8x10' 

33 21 4.2x10' 

50 39 2.3x10' 

150 20 4.4x10' 
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Table 6: Evolution of Surface Abundances (15 M©) 







H-Burn 


a 


He-Burn'' 


PreSN 






H-Burn 


a 


He-Burn^ 


PrcSN 




no- rot 


He 


U.ZlbZ 




0.2781 




U.o292 




12n 


/.OX iU 


-3 


2.3x10- 


-3 


i.6x iU 


-3 


rot 


^Hc 


0.2777 




n 0870 




0.3401 




12Q 


2.2x10- 


-3 


l.ZX lU 


-3 


9.1x10" 


-4 


rot+B 


^Hc 


0.2770 




n OSQQ 

U.zooo 




0.3365 




12c 


2.2x10^ 


-3 


1 /I v 1 ri" 

1.4X iU 


-3 


1.0x10" 


-3 


no- rot 


13 n 


Q /I X/ 1 n- 
0.4X iU 


-5 


1.3x10" 


-4 


i.UX lU 


-4 


15 m 

In 


Q o X/ 1 n- 
o.Z X lU 


-6 


1.8x10" 


-6 


l.o X iU 


-6 


i Ul 


13 p 




-4 


1.7x10" 


-4 


1 zLv 1 0" 


-1 






-6 


4.9x10- 


-7 




-7 


rotH-B 


13C 


1.3x10" 


-4 


1.9x10" 


-4 


1.5x10" 


-4 




1.6x10- 


-6 


4.9x10" 


-7 


4.3x10" 


-7 


no-rot 


14N 


8.2x10- 


-4 


1.3x10" 


-3 


3.3x10" 


-3 


160 


7.6x10- 


-3 


7.6x10- 


-3 


6.3x10- 


-3 


rot 


14N 


1.6x10' 


-3 


3.0x10" 


-3 


4.4x10" 


-3 


160 


7.5x10- 


-3 


7.0x10- 


~3 


6.0x10" 


-3 


rot+B 


14N 


1.5x10" 


-3 


2.6x10" 


-3 


4.1x10" 


-3 


160 


7.5x10" 


-3 


7.3x10" 


-3 


6.1x10" 


-3 


no-rot 


1^0 


3.1x10- 


-6 


3.6x10" 


-6 


5.7x10" 


-6 


180 


1.7x10- 


-5 


1.6x10- 


-5 


5.7x10- 


-6 


rot 




5.4x10" 


-6 


8.2x10" 


-6 


8.1x10" 


-6 


180 


1.5x10- 


-5 


1.0x10- 


-5 


7.6x10" 


-6 


rot-FB 


1^0 


4.9x10" 


-6 


7.0x10" 


-6 


7.2x10" 


-6 


180 


1.5x10- 


-5 


1.1x10- 


-5 


8.6x10" 


-6 


no-rot 


23Na 


3.4x10- 


-5 


3.6x10" 


-5 


6.4x10" 


-5 


19p 


4.1x10- 


-7 


4.0x10- 


-7 


3.2x10" 


-7 


rot 


23Na 


4.5x10" 


-5 


6.1x10" 


-5 


8.0x10" 


-5 


19p 


3.7x10- 


-7 


3.3x10- 


-7 


2.6x10" 


-7 


rot-FB 


23Na 


4.2x10" 


^5 


5.2x10" 


-5 


7.5x10" 


-5 


19p 


3.8x10- 


-7 


3.5x10- 


-7 


2.8x10" 


-7 


no-rot 


iiB 


3.8x10" 


-9 


2.6x10- 


10 


1.8x10- 


10 


9Bc 


1.7x10- 


10 


3.5x10" 


13 


2.4x10- 


13 


rot 


iiB 


2.4x10' 


10 


1.3x10- 


11 


8.4x10- 


12 


9Bc 


9.8x10- 


14 


6.0x10" 


16 


4.0x10- 


16 


rot-FB 


i^B 


1.1x10- 


10 


6.5x10- 


12 


4.5x10- 


12 


9Be 


2.9x10- 


15 


5.7x10" 


18 


3.9x10- 


18 



NOTE: "35 % central hydrogen mass fraction; ^50% central helium mass fraction; 
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2 4 6 8 10 12 

enclosed mass (solar masses) 




4 6 8 10 12 

enclosed mass (solar masses) 




02468 10 12 - 02468 10 

enclosed mass (solar masses) enclosed mass (solar masses) 



Fig. 1. — Magnetic field structure and angular momentum distribution for the standard 
15 M0 model at at hydrogen depletion (upper left), helium ignition (upper right), hehum 
depletion (lower left), and carbon ignition (lower right). See Table 1 for the definitions of 
these times. The shaded regions are those portions of the star that are convective. In those 
regions large diffusion coefficients for angular momentum lead to nearly rigid rotation in all 
but the latest stages of the evolution. 
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2 4 6 8 10 12 

enclosed moss (solor mosses) 




2 4 6 8 10 12 

enclosed moss (solor mosses) 



Fig. 2. — Specific angular momentum distribution for tlie standard 15 Mq model with (upper 
panel) and without (lower panel) magnetic fields at different evolution stages. See Table 1 
for the definitions of these times. 
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- pre-SN 

I depletion 




5 10 15 

enclosed moss (solar mosses) 




Fig. 3. — Same as Fig. 2 but for 25 stars. 



